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ABSTRACT 

We present an asteroseismic analysis of the helium atmosphere white dwarf (a DBV) recently found in the 
field of view of the Kepler satellite. We analyze the 5-mode pulsation spectrum that was produced based on one 
month of high cadence Kepler data. The pulsational characteristics of the star and the asteroseismic analysis 
strongly suggest that the star is hotter (29200 K) than the 24900 K suggested by model fits to the low S/N 
survey spectrum of the object. This result has profound and exciting implications for tests of the Standard 
Model of particle physics. Hot DBVs are expected to lose over half of their energy through the emission of 
plasmon neutrinos. Continuous monitoring of the star with the Kepler satellite over the course of 3 to 5 years 
is not only very likely to yield more modes to help constrain the asteroseismic fits, but also allow us to obtain 
a rate of change of any stable mode and therefore measure the emission of plasmon neutrinos. 
Subject headings: Dense matter — Elementary particles — Stars: oscillations — Stars: variables: general — 
white dwarfs 



1. INTRODUCTION 

Much to everyone's surprise considering the statistical 
odds, a helium atmosphere pulsating white dwarf (a DBV) 
was discovered in the 105 square degree part of the sky mon- 
itored by the Kepler satellite. GALEXJ192904.6+444708 
(a.k.a WD J 1929+4447 or KIC 8626021 in the Kepler input 
catalog) was found to be a good DBV candidate in a small 
auxiliary survey undertaken by Oestensen et al. (201 la), fol- 
lowing the null-detection of pulsations in the 17 WDs ob- 
served during the survey phase of the Kepler mission (Oes- 
tensen et al. 2010, 2011b). A low S/N survey spectrum 
of the object revealed its physical parameters to be Tcff = 
24900 ± 750 K and log g = 7.91 ± 0.07, placing it right in 
the middle of the DBV instability strip. One month of short 
cadence (58.8 s exposures) data was collected by the Kepler 
spacecraft in October and November 2010, revealing a pul- 
sation spectrum with five roughly equally spaced pulsation 
modes, three of which appears as evenly split triplets. The 
Fourier Transform and pulsational properties are presented in 
the discovery paper (Ostensen et al. 201 la). In this paper, we 
present the first asteroseismic analysis of WD J1929+4447. 

Our asteroseismic analysis strongly suggests that WD 
J 1929+4447 is more likely a blue edge DBV, with a pulsa- 
tion spectrum reminiscent of EC 20058-5234's (a.k.a QU Tel, 
Koen et al. 1995). This makes the star even more interesting 
as it then becomes the second known hot DBV with stable 
periods. Stable hot DBVs are good candidates to study the 
emission of plasmon neutrinos (Winget et al. 2004). An ob- 
served rate of change of the periods of any phase-stable mode, 
combined with models, will allow us to place constraints on 
the emission of plasmon neutrinos. This concept was at- 
tempted with EC 20058-5234 (Dalessio et al. 2010; Bischoff- 
Kim 2008), but without much success because of crowding 
of the field ai'ound the star. WD J 1929+4447 does not suf- 
fer from such pollution and may yield more positive results. 
Evolutionary rates of period changes have succesfully been 
measured for other pulsating white dwarfs, the most notable 



example being G117-B15A (e.g. Kepler et al. 2005). Hot 
DBVs cool a factor of a hundred times faster than a cooler 
white dwai-f Hke G117-B15A (Corsico & Althaus 2004). If 
it were not for the pollution from a field star for EC 20058- 
5234, an evolutionary rate of period change could have been 
obtained for that object based on a 5 year baseline of data. 

In section |2] we summarize the observed pulsational prop- 
erties of WD J1929+4447. In section [3741 we use these prop- 
erties to perform an asteroseismic study of that star and find 
a temperature of 29200 K, at odds with the spectroscopy. We 
discuss this discrepancy and other results in section |4] We 
summarize and conclude in section |5] 

2. PULSATIONAL PROPERTIES 

A total of 5 independent modes were found in the high 
signal-to-noise Fourier Transform (FT) of the data collected 
by the Kepler satellite. The 5 modes are listed in table [T]). 
We reproduced and preserve d the labels given to ea ch mode 
(first column in table [Tl from lOstensen et al]|201 lal) for con- 
sistency. The modes labeled /i o, /2,0i and /a q are the central 
members of clearly defined triplets in the FT. This allows us 
to positively identify these modes as £ = 1, m = modes 
and helps us constrain the asteroseismic fits. £ = 1 triplets are 
caused by rotat ion of the star. From the frequency splitting 
of these triplets. lOstensen et al J determine a rotation period of 
1.2 day for the star, typical for white dwarfs. Table [1] also 
shows the periods for our best fit model, discussed in section 

m 

As noted in the discovery paper, WD J 1929+4447 is pul- 
sationally most similar to EC 20058-5234, a blue edge DBV. 
Both stars have short period modes observed in their pulsa- 
tion spectrum. EC 20058-5234 has at least 6 independent 
modes under 400 seconds and 2 around 500 seconds (Sulli- 
van et al. 2008). But EC 20058-5234 has a temperature spec- 
troscopically and asteroseismically determined to be around 
28000 K (Beauchamp et al. 1999; Bischoff-Kim & Metcalfe 
2011), while Ostensen et al. (2011a) determined a tempera- 
ture of 24900 K for WD J1929+4447, placing it in the middle 
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TABLE 1 

Observed periods and the corresponding periods for the 

BEST fit model 



Mode 


Observed 


Amplitude 


Model 


t 




Label 


Period (s) 


(mma) 


Period (s) 






/2,0 


197.11 


2.92 


197.28 


1" 


3 


/l,0 


232.02 


5.21 


232.46 


1" 


4 


/3,0 


271.60 


1.87 


271.27 


1" 


5 


/4 


303.56 


1.29 


303.64 


1 


6 


/5 


376.10 


1.05 


375.87 


2 


15 



"These modes are identified as £ = 1, m = based on being the central mem- 
bers of well-defined triplets in the FT. 

of the DBV instability strip and not on the blue edge. 

3. ASTEROSEISMIC ANALYSIS 
3.1. The Models 

To compute all our models, we used the White Dwarf 
Evolution Code (WDEC). The WDEC evolves hot polytrope 
models from temperatures close to 100,000K down to the 
temperature of our choice. Models in the temperature range of 
interest for the present study are thermally relaxed solutions 
to the stellar structure equations. Each model we compute for 
our grids is the result of such an evolutionary sequence. 

The WDEC is described in detail in Lamb & van Horn 
(1975) and Wood (1990). We used smoother core composi- 
tion profiles and experimented with the more complex pro- 
files that result from stellar evolution calculations (Salaris et 
al. 1997). We updated the envelope equation of state tables 
from those calculated by Fontaine et al. (1977) to those given 
by Saumon et al. (1995). We use OPAL opacities (Iglesias & 
Rogers 1996) and plasmon neutrino rates published by Itoh et 
al. (1996). 

DBVs are younger than their cooler cousins the DAVs. 
Time dependent diffusion calculations (e.g. Dehner & 
Kawaler 1995) show that at 24000 K, a typical temperature 
for a DBV, the carbon has not fully settled into the core of 
the star yet. We expect the helium layer to be separated into 
a mixed He/C layer with a pure He layer on top, as shown in 
Fig. [U Following Metcalfe (2005), we adopted and parame- 
terized this structure in our models. 

We calculated grids of models using the WDEC and then 
ran a fitting subroutine to match the periods of the models 
(Peak) with the observed periods and calculated residuals us- 
ing the usual formula 



ctrms — 



Hobs 



(1) 



where nobs is the number of periods present in the pulsation 
spectrum. 

3.2. Parameters and grids 

In our asteroseismic fits, we vary up to 6 parameters; the 
effective temperature, the mass and 4 structure parameters. 
There are two parameters associated with the shapes of the 
oxygen (and carbon) core composition profiles: the central 
oxygen abundance (Xq), and the edge of the homogeneous 
carbon and oxygen core (q&n)- For envelope structure, Mcnv 
marks the location of the base of the helium layer and Mhc 
the location where the helium abundance rises to 1 (see Fig. 
[T]i- Mcnv and Mhc are mass coordinates, defined as e.g. 
Mcnv = -log(l - M(r)/M,), where M(r) is the mass en- 
closed in radius r and is the stellar mass. 
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Fig. 1. — Helium mass fraction of model DBs. The symbols show the 
profile that results from the time dependent diffusion calculations of Dehner 
& Kawaler (1995). The solid line shows the profile for our best asteroseismic 
fit of WD J1929+4447 (see section[34l. 

With only 5 periods, 6 parameter fits are underconstrained. 
While there is much to learn from attempting 6 parameter fits, 
we settled for 4 parameters (later relaxing a 5*''). There are 
three logical ways of fixing 2 of the 6 parameters: 1) fix the 
mass and effective temperature according to the spectroscopy, 
2) fix the core structure parameters, and 3) fix the envelope 
parameters. In light of the basic pulsational properties of the 
object (see also section [33] l. option 1) seems the least promis- 
ing. Of the other two choices, we decided to first try fixing 
the envelope structure to that of Dehner & Kawaler (1995). 
This meant that Mcnv was fixed to -2.80 for all models and 
Mhc was first fixed to -5.50 for models between 20000 K and 
22000 K, -5.90 for models between 22200 K and 26000 K 
and -6.30 for models between 26200 K and 28000 K, to match 
the profiles resulting from the diffusive settling computed by 
Dehner & Kawaler (1995). 

We calculated 2 grids of models; one covering a broad 
range of parameter space at lower resolution and one zoom- 
ing in on the promising area of parameter space with higher 
resolution. The range of parameters and step sizes for each 
are summarized in table |2] The broad grid contains 305,820 
models that succesfully converged (out of 315,700 total com- 
puted) and the higher resolution grid contains 403,065 con- 
verged models, with a similar success rate. Table|2]also shows 
the parameters for the best fit models, discussed in section [34l 

3.3. Average period spacing 

The average period spacing provides a useful asteroseismic 
measure of the mass and temperature of the star, indepen- 
dent of the details of internal chemical composition profiles. 
Higher k modes are not trapped in the core and according 
to asymptotic theory, they are evenly spaced in period. This 
spacing is given by (Unno et al. 1989) 



AP: 



+ 1)1/2 



ri 



'■2 N , 
— dr 
r 



1 -1 



(2) 



where ri and r2 are turning points of the mode and N is the 
Brunt- Vaisala frequency. The asymptotic period spacing is £ 
dependent, with higher £ modes having smaller spacing. The 
dependence on the Brunt-Vaisala frequency leads to higher 
mass and higher temperature models having a smaller period 
spacing. 

In the case of WD J1929H-4447, we must be careful when 
using average period spacing arguments. First, the modes 
are not higher k modes, so at least a subset of them may 
be strongly trapped in the core (we are not in the asymptotic 
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TABLE 2 

Parameters covered in the model grids and best fit parameters for WD J 1929+4447 



Parameter 


Grid 1 

Minimum Maximum Step Size 


Minimum 


Grid 2 
Maximum 


Step Size 


Best fit 


Toft (K) 


20000 28000 200 


24100 


30000 


100 


29200 


M, (Solar) 


0.500 0.690 0.010 


0.555 


0.650 


0.05 


0.570 


Menv 


-2.80 (fixed) 




-2.80 (fixed) 




-2.80 


Mho 


-5.50, -5.90, or -6.30 (see text) 


-6.10 


-6.50 


0.20 


-6.30 


Xo 


0.00 1.00 0.10 


0.40 


1.00 


0.05 


0.60, 0.65" 


qfm (M.) 


0.10 0.80 0.02 


0.32 


0.48 


0.02 


0.36 



"Close tie 



limit). Second, we only have a few modes to give us an av- 
erage period spacing. That being said and somewhat ignoring 
this warning, we proceed. 

From the first 4 modes in table [T] we determine an average 
period spacing of 35.9 seconds. This is to compare to 36.5 
seconds for EC 20058-5 234 (~ 28000 K) and 38.8 seconds 
for GD358 (- 25000 K, [Provencal et ai1l2009l) . The former 
marks the current blue edge of the DB V instability strip, while 
the latter is the prototype of a middle of the instability strip 
DBV. From these numbers alone we might conclude that WD 
J1929H-4447 is closer in temperature (and mass) to EC 20058- 
5234 than it is to GD358, at odds with the spectroscopy. 

To gain a better notion of how a 3 second difference in av- 
erage period spacing translates into a difference in mass and 
effective temperature, we calculated asymptotic period spac- 
ings for the models in our grids and pulled models with pe- 
riod spacings similar to that of EC 20058-5234 and another 
set with spacings similar to that of GD358 (performing some 
normalizing in resolution to merge grid 1 and grid 2 - 31500 
models total). We show the effective temperature and mass 
distribution of the two populations in Fig. |2] The 36.5 sec- 
ond period spacing population peaks at an effective tempera- 
ture of 27500 K while the other peaks at 25500 K, consistent 
with the effective temperatures of the two stars. This suggests 
that WD J1929H-4447, with its shorter period spacing, is hotter 
than 24900 K. 

3.4. Results of the astemseismic fitting 

The parameters of the two best fit models that resulted from 
our grid search are listed in table |2] and the periods of the 
marginally better one (Xq ~ 0.60) in table[T] For these mod- 
els, (Trms = 0.28 s. To get a sense of the quality of this 
fit, it is useful to compare it to similar fits performed in the 
litterature. The difficulty in comparing the results of the as- 
teroseismic analysis of different pulsating white dwarfs is that 
the number of observed periods to fit and the number of pa- 
rameters used varies. One way to compare them is to use the 
BIC parameter (Bayes Information Criterion). The BIC pa- 
rameter measures an absolute quality of the fit, by taking into 
account differing nurn bers of data points and free parameters 
(iKoen & Laneyll2000h . It penalizes fits involving a greater 
number of parameters relative to the number of data points. 
Applied to asteroseismology, it is given by 



BIC 



In(nobs) 

Hobs 



iiobs In(CT^Ms), 



(3) 



where Upar is the number of free parameters in the fit. A 
smaller BIC parameter indicates a better fit. 

For the present case, iiobs ~ 5 and iipai = 5. With 
CRMS — 0.28 s, this translates to BIC = —0.41. Among the 
best recent asteroseismic fits (Bischoff-Kim et al. 2008; Cas- 
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Fig. 2. — Temperature (upper panel) and mass (lower panel) distribution of 
a subset of model DBs. The thick red lined histograms show the distribution 
for models that have an asymptotic period spacing between 36 and 37 seconds 
(like EC 20058-5234) and the thin blue lined histograms show the distribution 
for asymptotic period spacings between 38.2 and 39.2 seconds (like GD358). 

tanheira & Kepler 2008), using the BIC, only one fit among 
6 comes out superior But that is a fit to a single mode star 
using 4 parameters. It is difficult to argue for the significance 
of such a fit. 

It is illuminating to also look at the parameters (mass and 
effective temperature) of models that did not fit quite as well, 
but were still good fits. Fig. [3] is a plot of where all the 
models that fit to better than (t-rms = 1-0 s lie in param- 
eter space. The best fit models follow a negatively sloping 
trend in the mass versus effective temperature plane. This is 
a direct consequence of Eq. [2] To fit a fixed set of observed 
periods, higher temperature models must compensate with a 
lower mass in order to keep the appropriate average period 
spacing. This trend frequently shows up in white dwarf aster- 
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Fig. 3. — Best fit models in 3 different planes of parameter space. The filled circles are models that fit the observed period spectrum with (TftMS 5: 1-0 s. The 
best fit models at 29200 K is circled. In the top panel (mass versus effective temperature), we also mark the spectroscopic value (diamond) and the formal 1 -a 
fitting error region around it. 



oseismology (e.g. Bischoff-Kim et al. 2008; Castanheira & 
Kepler 2008). The band of acceptable model fits in that plane 
misses the l-cr spectroscopic box completely. However, we 
do not consider this a failure as it is well known that the He 
1 lines in the DB instability region change only weakly with 
temperature. Thus, small calibration effects can easily throw 
the temperature by many times the formal 1-sigma fitting er- 
ror The mass of our best fit models, on the other hand, corre- 
sponds well with the mass determined from the spectroscopic 
surface gravity. 

Here it is worth summarizing our method and emphasizing 
one aspect of it before presenting one more result. We started 
out by fixing the helium envelope parameters to the values 



that fit the time dependent diffusion calculations of Dehner & 
Kawaler (1995) but when we zoomed in on the promising part 
of parameter space (grid 2 in table |2]i, we relaxed the param- 
eter that sets the thickness of the pure helium envelope. We 
also never fixed the effective temperature and treated it as a 
free parameter within a reasonable range. We now come back 
to Fig. [T] which actually shows quite an interesting result. 
On that graph, we show the helium abundance profile for our 
best fit model along with the profile calculated by Dehner & 
Kawaler (1995). Two things match beautifully: 1) The pro- 
files themselves (of course the base of the mixed He/C layer 
matches perfectly because we fixed that, but the base of the 
pure helium layer was eventually allowed to vary) and 2) The 
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effective temperature. The latter was allowed to vary between 
24100 K and 30000 K and settled on 29200 K. 

4. DISCUSSION 

This last result may be evidence that the calculations of 
time dependent diffusion calculations of Dehner & Kawaler 
(1995), as well as others since then e.g. Althaus et al. (2005) 
describe the evolution of PG1159 stars well. However, we 
also need to look at what the core structure is. Does it agree 
with stellar evolution calculations? The answer there is not 
quite as positive. We find a central oxygen abundance of 
0.60. With the standard reaction rates and treatment for con- 
vection, we expect the central abundance for a white dwarf of 
mass 0.550 M© to be between ^ 0.8 and ^ 0.9 (Althaus et 
al. 2010; Salaris et al. 1997). Perhaps more importantly, as 
our models are highly sensitive to that parameter (Bischoff- 
Kim & Metcalfe 2011), the location of the edge of the ho- 
mogeneous core is further inward than what stellar evolution 
calculations find. 

In this paper, we took three different approaches using 
the pulsations to determine an effective temperature for WD 
J1929H-4447: 1) Through a quick inspection of the pulsation 
spectrum, noting that low period modes were observed, 2) us- 
ing average period spacing arguments, and 3) by performing 
asteroseismic fits of the period spectrum. All three point to 
an effective temperature more similar to EC 20058-5234's 
and inconsistent with the spectroscopic value of Tch = 
24900 ± 750 K. A well-calibrated high-S/N spectrum is ur- 
gently needed to settle the temperature issue. 

5. CONCLUSIONS AND FUTURE WORK 

We performed an asteroseismic analysis of WD 
J1929H-4447, a DBV discovered in the field of view of 



the Kepler satellite. We find strong evidence that the star 
is a hot DBV. Our models also support the time dependent 
diffusion calculations of Dehner & Kawaler (1995), though 
more positive results from the study of other DBVs are 
needed to further support these results. Also, it would be 
worth approaching the problem from another angle. In 
section [32] we made the somewhat arbitrary choice of fixing 
the envelope parameters and allowing the core parameters 
to vary. We would most likely learn a lot by trying to fix 
the core parameters to those expected from stellar evolution 
(e.g. Althaus et al. 2010; Salaris et al. 1997) and allow the 
envelope parameters to freely vary over the entire reasonable 
range. Continuous observations with the Kepler satellite may 
also reveal more modes, hidden in noise at the current S/N 
level, allowing 6 parameter fits. 

The discovery of a hot DBV is significant, as it gives us an- 
other chance to study plasmon neutrino emission from white 
dwarfs. At 29200 K, WD J1929H-4447 should be radiating as 
much energy through the emission of photons as through the 
emission of plasmon neutrinos (Winget et al. 2004). The extra 
cooling due to neutrinos is therefore easily detectable if one 
could measure the cooling rate for the star with any signifi- 
cance. An observed period change in any of the modes can 
be a reliable estimator for the cooling rate, and only requires 
that at least one of the observed modes are phase-stable over a 
sufficient number of years to detect the variation. If Kepler is 
permitted to maintain observations of WD J1929H-4447 for a 
period of five years, the cooling rate should be determined to 
a sufficient precision to establish the respective contributions 
from thermal radiation and plasmon neutrinos. 
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